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ABSTRACT 

The major mass fraction of the envelope of hot luminous stars is radiatively stable. However, the partial 
ionisation of hydrogen, helium and iron gives rise to extended sub-surface convection zones in all of them. In 
this work, we investigate the effect of the pressure induced by the turbulent motion in these zones based on 
the mixing length theory, and search for observable consequences. We find that the turbulent pressure fraction 
can amount up to ~ 5% in OB supergiants, and to ~ 30% in cooler supergiants. The resulting structural 
changes are, however, not significantly affecting the evolutionary tracks compared to previous calculations. 
Instead, a comparison of macroturbulent velocities derived from high quality spectra of OB stars with the 
turbulent pressure fraction obtained in corresponding stellar models reveals a strong correlation of these two 
quantities. We discuss a possible physical connection, and conclude that turbulent pressure fluctuations may 
drive high-order oscillations, which — as conjectured earlier — manifest themselves as macroturbulence in the 
photospheres of hot luminous stars. 

Subject headings: stars: massive, stars: evolution, turbulence, convection, stars: oscillation, line: profile 


1. INTRODUCTION 


the interstellar medium (Kobayashi et al. 2006 

i, for regulat- 
, and as pro- 

ing star formation (Mac Low & Kies sen 2004 

genitors of supernovae and gamma ray bursts ( 

^angerf2012). 


However, various physical processes acting in massive stars 
are as yet not well understood, therefore preventing the cur¬ 
rent stellar models to self-consistently explain several general 
observational properties, such as the mass-discrepancy (Her- 
rero et ak|l992| Markova & Puls||2015)), the presence of the 


Humphreys-Davidson (HD) limit (Humphreys & Davidson 

7979}, and the position of the term inal-age- main-sequence 
in the HR diagram (e.g., |Vink et ~aT7| |2010| |Castro et al.| 
2014}. Radiation pressure dominated layers, envelope infla- 
tion, clumpy stellar winds, and dynamical instabilities are fea¬ 
tures which commonly occur in these objects, but which are 
not yet thoroughly investigated. Some of these phenomena 
have been connected to the presence of turbulent motions in 
the outer layers (e.g., |de Jager|[ 1984} |Cantiello et al.l|2009} 
Despite its s implicity, the mixing-length theory (MLT, Bohm- 
Vitens ep 958) 1 for convection has been very successful in de¬ 
scribing the main features of non-adiabatic turbulent stellar 
envelopes (e.g., |Trampedach et al.|2014} . However, the MLT 
is expected to have shortcomings when the convective veloc¬ 
ities approach the speed of sound , as it is expected in the en¬ 
velopes of very luminous stars (Sanyal et al. 2015) . Inparticu- 
lar, it has been argued by de Jager (1984) and Maeder) (12009 1 
that turbulent pressure may become important in this situa¬ 
tion. It is the aim of the present paper to discuss the role of 
turbulent pressure in the envelopes of massive stars, investi¬ 
gate its effects on the stellar stmcture, and focus on its poten¬ 
tial observational signatures. 

2. METHOD 

We use the Lagrangian one-dimensional hydrodynamic 
stellar evolution code BEC (Heger et al.|2000[ Petrovic et al. 


2005 1 Brott et ak|201 1 ) for computing massive star models. 
It treats convection following the MLT with a mixing-length 
parameter of a — 1.5 (Bro tt et al.|201 1 1 . The opacity of stellar 
matter is interpo lated from the OPAL opacity tables (Iglesias 
& Rogers|1996 ), and mas s-loss by stellar wi nd is adopted fol- 


lowing the prescription of|Vink et al.| 

We modified the BEC code in order to include the turbu¬ 
lent pressure term in the stellar envelopes, in particular in 
the convective zones associated with the iron opacity peak 
(FeCZ), and with the partial ionisation of helium (HeCZ), 
and hydrogen (HCZ). We neglect the turbulent pressure in the 
convective core, where it is typically more than 7 orders of 
magnitude smaller than the ideal gas pressure. We included 
the turbulent pressure P tul -b and the turbulent energy density 
e tul b in the momentum, energy transport and energy conserva- 


1991 1 [Jiang & Huang 

1997 

Trampedach et al.|20l4 

i: 


1997 Stothers 2003 Maeder 2009 


Pturb = Zpvj. 


_ 3 Pturb 

Cturb ~ ~ -> 

2 p 


( 1 ) 


where p is the local mass density, if i s a p arameter chosen to 
be if = 1/3 for isotropic turbulence ( |Stothers||2003) |Maeder| 
2009), and \y is the local convective velocity. As the convec¬ 
tive velocities approach the local sound speed, we limit the 
convective velocities to this value, i.e. 


v c < c s 



pm H 


p 

1 gas 

P 


( 2 ) 


where T is the local temperature, kb is the Boltzmann con¬ 
stant, p is the mean molecular weight, and m H is the proton 
mass. Here, we use the isothermal sound speed, as the lay¬ 
ers under consideration are characterised by a small ratio of 
the local thermal-to-dynamical time scale (r,;, lTd yn < 1)- As 
a consequence of Eq[2] the turbulent pressure in our models 
can not exceed a value of one third of the local gas pressure. 
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Fig. 1.— Stellar evolution tracks calculated including turbulent pressure and 
turbulent energy density, indicated by colored dots. Superposed are the tracks 
by|Brott et al.|(20l 1 \ (black lines). The colour indicates the maximum frac¬ 
tion of turbulent pressure occurring within each stellar model (see color bar 
to the right), where fractions above 3% are indicated in red. Stellar models at 
Teff < 10 kK coloured in red show contributions from P tur b up to ss 33% and 
hot luminous OB stars coloured red up to a 5%. The stellar masses next to 
the tracks indicate the initial mass of the models. The da shed line indicates 
the position of the HD-limit (Humphreys & Davidson] 1 9791 . 

3. RESULTS 

We computed a set of stellar models and stellar evolutionary 
tracks for non-rotating stars with Galactic metallicity in the 
initial mass range 7-80M G , including the physics described 
in Sect . 2. All input and physics parameters were chosen as 
in |Brott et al.| ( |20lT) . 


3.1. Evolutionary Tracks 

The evolutionary tracks of our newly computed models are 
plotted in Fig. [I] and superposed to those obtained by Brott 


et al. (20111. Throughout most of the evolution, the effect 


of turbulent pressure and energy on the luminosity and the 
surface temperature of the models is small, such that our 


new evolutionary tracks nearly coincide with those from Brott 
et al. ( |2011| >. According to the colour scheme in Fig. [T| 


which indicates the maximum fraction of turbulent pressure 
inside the stellar models, the HR-diagram can be roughly di¬ 
vided into three areas: the hot (T e ff> 10 4 K) and luminous 
(log(L/Lo) > 4.5) stars, the cool stars (T e ff < 10 4 K) of any 
luminosity, and the hot stars with log(L/Lo) < 4.5. 

In the latter regime, the turbulent pressure does not exceed 
a fraction of the total pressure of a few tenths of a percent 
anywhere inside the stellar models. On one hand, the iron 
opacity peak in these models is located deep inside the enve¬ 
lope in nearly adiabatic layers, which leads to relatively low 
convection velocities. On the other hand, the models are too 
hot to contain a hydrogen ionisation zone. Finally, they do 
show a Hell partial ionisation zone close to the surface, which 
is, however, not vigorous enough to play a significant role 
( |Cantiello et al.[20091 1. 

In the cool supergiant region, the turbulent pressure fraction 


inside the star can be very significant. The maximum value is 
typically 25-30%, while the highest possible value, 33%, is 
reached for the most luminous models. The computed tran¬ 
sonic convective velocities, of the order of 20-30 km/s, arise 
in the HCZ very close to the surface, where high opacities 
induce high local Eddington factors, density inversions, and 
high degrees of superadiabaticity (Sanyal et al .|[2015) . Such a 
high turbulent pressure fraction is possible here since ideal gas 
pressure is by far the dominant contribution to the total pres¬ 
sure in the envelopes of the cool supergiants. Even though a 
significant fraction of the total pressure arises from the turbu¬ 
lent motion, the evolutionary tracks are only shifted by few 
tens of Kelvin towards lower temperatures, compared to the 
tracks without turbulent pressure, which corresponds to a ra¬ 
dius increase by a few percent. This is due to the fact that the 
region which contains high turbulent pressure contains very 
little mass (see below). 

In the hot and very luminous stars, the turbulent pressure 
can account for up to ~5% of the total pressure. The evo¬ 
lutionary tracks are not significantly different from those of 
Brott et al. (2011), showing displacements of the order of hun¬ 


dreds of Kelvin towards lower effective temperatures. Tur¬ 
bulent pressure becomes more important as the stars expand 
during their main sequence evolution, as the convective ve¬ 
locities increase. It achieves a maximum in the O super¬ 
giant regime, and then decreases in the B supergiant regime 
as the iron opacity peak moves deeper inside, but turbulent 
pressure remains significant for surface temperature above 
~10000K. The maximum turbulent pressure in these models 
occurs within the FeCZ, where high local Eddington-factors 
are achieved, giving rise to envelope inflation and density in¬ 
versions (Sanyal et al. [2015) . It does not reach as high frac¬ 
tions as in the cooFsupergiant because of the predominance 
of radiation pressure in the hot star envelopes. 


3.2. The structure of the outer layers 

Following the results shown in Sec 3.1, we investigate 
the structure of the 20 M Q and 60 M Q models in more detail. 
Figs. [2] and [3] show the relative fraction of the turbulent pres¬ 
sure as function of the optical depth in the stellar envelopes of 
both models throughout the evolution. 

In the 20 M 0 model (Fig. [ 2 }, we find that for T e g> 10 4 K 
the turbulent pressure fraction has its maximum in the FeCZ, 
located at an optical depth of log(r)« 2.5 - 4. The turbu¬ 
lent pressure fraction is about 1.8% at 7 d r - 25 000K, and 
the FeCZ moves deeper inside the star as the stellar model 
expands during its evolution. Once the stellar model reaches 
effective temperatures well below 10000K, the HCZ arises. 
This convective region, within which the maximum turbulent 
pressure fraction arrives at about 25-30%, reaches the stellar 
surface in stellar models with T e ff < 8000 K, and covers an 
extended range of optical depths. The bottom panel of Fig. [2] 
compares the density profiles and the radial extent of the con¬ 
vection zones for a model with 7 d r — 7000 K and turbulent 
pressure included to the same model where turbulent pressure 
was disregarded. 

In the envelope of the 60 M G model, we find convective 
zones associated to the partial ionisation of iron and helium, 
while the model does not become cool enough to show hydro¬ 
gen recombination. 

As shown in Fig. [3] the FeCZ is located at an optical 
depth of log(r) > 1.5Twith P tur b/P rising to ~5% for 7' d r - 
40 000 K. As the star evolves the convective region moves 
deeper inside the star. Figure [3] shows that due to the tur- 
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Fig. 2.— Top: ratio of turbulent-to-total pressure (color coded) as a func¬ 
tion of effective temperature and optical depth throughout the evolution of 
the 20 M 0 model. The FeCZ is visible in the upper-right comer, the HeCZ 
lies across the diagram, and the HCZ appeals for r e ff below ~ 10 000 K. The 
dashed line indicates the selected density profile shown in the bottom panel. 
Bottom: Comparison of the density (in gem -3 ) as a function of the radial 
coordinate R in the outer layers of our 20 M© star at T^ff - 7000 K (blue 
solid line) with that in an identical model where the turbulent pressure term 
was switched off (black dashed line). The shaded dark-red regions confined 
within the dotted vertical lines denote the HeCZ and HCZ in the model cal¬ 
culated without P(urb ? while the light-red regions correspond to the same con¬ 
vection zones in the model calculated with P tur b- The numbers associated to 
the arrows on the top of the figure indicate the optical depth log(r) of the bor¬ 
ders of the convective zones in the case with P tur bincluded (with S indicating 
the surface). 

bulent pressure, the star increases its radius by few percent, 
which leads to a slightly reduced density and an increase of 
the radial extent of the iron convection zone. 


4. COMPARISON WITH OBSERVATIONS 

In particular the iron opacity peak at T »20 0000 K is 
known to induce a variety of dynamical phenomena at the 
stellar surface, e. g. pulsations in /j-Cephei and slowly pulsat- 
ing B stars (SPB, jMiglio et al.|2007j [Pamyatnykh] 1999}), and 
stochastically excited travelling waves generated by turbulent 
motions in the FeCZ (|Goldreich & Kumar||l990t Cantiello 


et al. 

2009 

Belkacem et al.|2010; Samadi et al. (2010[ [Mathis 

et al. 

2014 

) leading to a subsonic small scale velocity held at 


Fig. 3.— Same as Fig. [2] but for the 60 M© model. Top: ratio of turbulent-to- 
total pressure (color coded) as a function of effective temperature and optical 
depth during the evolution of the 60 M© model, showing the FeCZ and HeCZ 
are at higher and lower optical depths, respectively. 

Bottom: density profile and extent of the FeCZ at T e ff ss 40 000 K. 


the surface that has been proposed to be the physical origin of 
the so called “microturbulence”. 

The spectra of luminous O-B stars is known to be also 
affected by the so called macroturbulent broadening, an ex¬ 
tra line-broadening usually ad-hoc associated with large scale 
(compared to the line forming region) motions at the surface 
(see |Markova et al.| ( |2014| >, |Simon-Dfaz & Herrero| ( |2014| ), 
and references therein). Similarly to the case of microturbu¬ 
lence, convection might play a significant role in the origin of 
macroturbulence as well. This view is supported by the work 
of Sundqvist et al. (]2013|>, who showed that macroturbulence 
is suppressed in strongly magnetic massive stars, where the 
magnetic held is expected to at least partially inhibit convec¬ 
tion. A similar effect was found for intermediate mass mag¬ 
netic chemically peculiar stars ( |Ryabchikova et al.|1997) > and 
for spots in late-type stars (Strassmeier 2009). In this context, 
vigorous envelope convection in the temperature range of the 
hydrogen recombination (Fig. [2|» may be responsible for the 
non-thermal (macroturbulent) broadening observed in red su¬ 
pergiants ( Collet et al.|2007||Carney et al.|2008] l. 
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With this motivation, we pursue the hypothesis that the rel¬ 
ative strength of turbulent pressure in the sub-surface convec¬ 
tive zones is related to the appearance and strength of macro¬ 
turbulence at the stellar surface. We investigate the case of the 
luminous O-B stars, where the turbulent pressure constitutes 
up to ~5% of the total pressure in the FeCZ. 

We make use of the results from the quantitative spectro¬ 
scopic analysis of the rich s amp le of spec tra compiled by the 
IACOB project (Simon-Dfaz et al. |2()11 ||2015 1 . In particu¬ 
lar we benefit from the derived values of surface temperature 
(T e ff ), gravity (log(g)), projected rotational velocity (vsin(i)) 
and macroturbulent veloc ity (v mnrr n) fo r _a_ sample of ~ 300 


Galactic O-B stars used in Simon-Dfaz (2015 

In Fig. HI we compare the observationally derived macrotur¬ 
bulent velocities to the maximum fraction of turbulent pres- 
sure in our models in the spectroscopic HRD (sHR, |Langer| 
& Kudritzki||2014|. Stars presenting a clear signature of 
macroturbulence in their line profiles (i.e., a v macro larger than 
50km/s) are marked by a bigger black circle in Fig. [4] Inter¬ 
estingly, those stars are located mainly in regions of the sHR 
where the turbulent pressure yields the highest contribution to 
the total pressure in the FeCZ. At places where the models 
predict a small contribution from the turbulent pressure only 
very few stars show an unambiguously high macroturbulent 
velocity. We discuss these exceptions at low ££ in the next 
section. 

The agreement becomes even more striking looking at 
Fig. [5] where the observed macroturbulent velocities are di¬ 
rectly plotted against the maximum fraction of turbulent pres¬ 
sure in our models. This plot reveals a clear correlation be¬ 
tween macroturbulent velocity and turbulent pressure starting 
from /P ~ 0.005, with a Spearman’s rank correlation co¬ 
efficient of 0.812. 

5. THE CONNECTION TO HIGH-ORDER PULSATIONS 

Figure[4]revealed about ten stars with log(_S?/»2o) < 3.3 for 
which our models predict very small turbulent pressure con¬ 
tributions, but which unambiguously show macroturbulence 
at a level above 50km/s (see also Fig. [5|. Interestingly, we 
find that all these stars are located inside or very near to the 
region where stars are expected to be pulsationally unstable 
to high-order g-modes (|Miglio et al.|2007] Pamyatnykh| 1999 


see also Simon-Dfaz 2015|), which is drawn as a gray band in 


Fig- 0 Indeed, |Aerts et ah] ( |2009) > showed that the collective 


effect of high-order non-radial pulsations may produce a ve¬ 
locity field in the spectra of hot stars which resembles closely 
what is observed as macroturbulence in the stars discussed in 
Sect. 4. 

Consequently, we interpret these stars as affected by a 
macroturbulent broadening that can be explained via a heat- 
driven pulsational origin. However, the homogeneity of the 
spectroscopic signatures of macroturbulence over the whole 
effective temperature range in the luminous stars calls for a 
single dominant mechanism to produce it (Simon-Dfaz et al. 
in prep.). This is not the case when considering only classical 
(/(-mechanism) instability domains, which do not cover the 
full region where most of the stars showing a macroturbulent 


velocity field are observed (see 

Simon-Dfaz 2015 Godart et 

al. in prep.). Shiode et al. (2013 

) and Aerts & Rogers (2015) 


consider gravity-waves originating in the convective core as 

1 The results of the analysis of a much larger dataset, along with a thorough 
empirical description of the behaviour of macroturbulent broadening in the 
whole O-B star domain will be presented in Simon-Dfaz et al. (in prep). 


the cause of macroturbulence. Whereas |Shiode et al.| ( |2013] > 
based on their massive star models find that the surface ve¬ 
locity fluctuatio ns do not exceed 10 m/ s even in their most 
optimistic case, Aerts & Rogers ( |2015) >, based on 2-D non¬ 
linear simulations of convection-driven waves in a modified 
3 M 0 model concluded that it might explain the macroturbu¬ 
lence observed in O-stars. 

Within our scenario, the Reynolds stresses associated with 
turbulent pressure induce uncorrelated turbulent pressure fluc¬ 
tuations in the form 


SPturb ~ pv" c , 


(3) 


( Goldreich & Kumarj 1990j Grigahcene et al.|2005 Lecoanet 


|& Quataert|201 3[ |Shiode et al.|2013| ) where 5P tur b is the La- 
grangian pressure perturbation associated with the convec¬ 
tive motions. Such stochastic fluctuations at the percent level 
can produce a strong local deviation from hydrostatic equilib¬ 
rium, and will thus naturally excite high-order pulsations in 
the range of eigenmodes, which are closest to the spectrum of 
the fluctuations. 

This would imply that also in the luminous stars with 
log(^/if 0 ) > 3.3 in Fig. J4] the macroturbulence may be 
signifying high-order pulsations, which are, however, excited 
by turbulent pressure fluctuations rather than through the k- 
mechanism or by strange mode instability. If so, we can on 
one hand expect that linear pulsation analyses which include 
the R ey nolds stress tensor, as e.g., in Dupret et al. (|2004 


2005) or|Antoci et al.|(|2014|l, may uncover that stars in a large 
fraction of the red and orange coloured region in Fig. [4] are 
unstable to high-order g-mode pulsations. On the other hand, 
as the pressure fluctuations in these stars, as predicted by our 
simple analysis, can amount up to 5% of the total pressure, it 
is conceivable that in linear stability analyses, which require 
the growth of the instability from infinitesimal perturbations, 
an instability is not detected in all stellar models where high- 
order g-mode can be excited through finite amplitude pressure 
perturbations. 

6. CONCLUSIONS 

We implemented the effect of the turbulent motion of con¬ 
vective eddies in a simple formalism in the momentum and 
energy equations of our stellar evolution code. By comparing 
to previous computations (Brott et al. |201 1 1, we find that the 
turbulent pressure does not alter the stellar structure signifi¬ 
cantly. However, we find maximum turbulent pressure con¬ 
tribution of up to 5% and 30% in our models for O-B super¬ 
giants and cool red supergiants, respectively. By comparing 
the maximum turbulent pressure contribution in our models 
with spectroscopically derived macroturbulent velocities for 
a large sample of Galactic OB stars (Sim6n-D faz]|2015] l, we 
find both quantities to be strongly correlated. 

Several less luminous stars, in which the turbulent pres¬ 
sure is thought to be small, show nevertheless high macro¬ 
turbulent velocities. These are located close to or inside the 
region where linear pulsation analysis predicts high-order g- 
mode pulsations, arguing therefore for /(--mechanism pulsa¬ 
tions, and not turbulent pressure fluctuations, as the origin of 
the macrot urbulence phenomenon, in line wi th previous sug¬ 
gestions ( |Aerts et al.|2009 Sim6n-Dfaz|2015| l. 

We argue that the turbulent pressure fluctuations in hot 
luminous stars can excite such high-order pulsations, most 
likely non-radial g-modes, which may explain the occurrence 
of macroturbulence in stars which are found outside of the 
currently predicted pulsational instability domains. This view 
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Fig. 4. — Spectroscopic HRD with coloured regions representing the maximum fraction of tu rbulent pressure as derived from a best-fit of the tracks in Fig. IT] 
Circles represent the observed O-B stars, colour coded according to their spectral line shape as in|Sim6n-Dfaz|{20151. i.e. following the ratio between the projected 
rotational velocity v sin i and the macroturbulence velocity v macro ■ In the left panel stars with line profiles showing a clear contribution from macroturbulence 
are located, while the right panel includes stars with line profile showing a dominant rotational broadening and a less clear contribution from macroturbulence. 
The bigger circles bordered in black indicate stars showing v macT „ higher than 50 km/s. The gray bands indicate the K-mechanism instability strip for SPB stars 
(Might) et al.|2007) . 



Fig. 5.— v macro as a function of the maximum fraction of turbulent-to-total 
pressure derived from a fitting of the tracks described in Sect. 3. The circles 
are color coded as in Fig.Rlwith the crossed-black circled dots indicating the 
stars with v„ mcro > 50 km/s and log< 3.3. The dashed line indicates 
the 50 km/s level (see Fig.QJ. 
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is in agreement with the indication that macroturbulence can 
be suppressed in strongly magnetic stars, given that such a 
field might effectively inhibit convective motions in the FeCZ 
(Sundqvist et al. 2013) l. 

At the moment, turbulent pressure fluctuations appear to 
be the only mechanism which may excite high-order oscil¬ 
lations in luminous stars (log(L/Lo) > 4.5) in the wide effec¬ 
tive temperature regime for which strong macroturbulence is 
observed. 
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